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Abstract 


The protection of metal parts from corrosion is important to save the cost 
of production of goods and lives. Twelve triazolopyrimidinone 
derivatives were investigated for their ability to inhibit the corrosion of 
metal via density functional theory and Monte Carlo (MC) simulation 
approaches. The electronic properties and global and local reactivity 
descriptors were computed at the B3LYP/6-31G(d) level of theory, using 
Spartan 14 software. The Mulliken charge distribution and electrostatic 
potential (ESP) maps were used to locate the reactive sites. The adsorption 
of the compounds on the surface of the stable Fe(110) was investigated in 
an aqueous solution using the adsorption locator module on Material 
Studio software. The values of energies of their frontier molecular orbitals 
and other reactivity descriptors suggest that these molecules displayed 
good corrosion inhibition potentials. The Mulliken charge distribution and 
electrostatic potential maps showed that they can donate electrons to the 
metal and at the same time accept electrons via back-donation, also 
rationalized by the values of their electron back-donation calculated. The 
molecules interacted well with Fe(110) surface in an aqueous medium, as 
revealed by the MC results. 


y) © This work is licensed under the Creative Commons Attribution-Non- 
wy Commercial 4.0 International License. 


Introduction 


The use of metals in construction materials, oil and 
gas, and petroleum industries is very common. 
They, however, readily undergo degradation via 
chemical or electrochemical reactions with their 
surroundings [1], this is because they are 
thermodynamically unstable in their pure forms. 
Corrosion is a major issue that has affected several 
industries worldwide, destabilizing the world 
economy. It has resulted in the breakage of oil and 
underground water pipes, failure of bridges, 
contamination of products from leakages of pipes, 
damage of metallic equipment, and contamination 
of air through the release of harmful pollutants 
from iron corrosion [2]. In 2013, the cost of 
corrosion globally was estimated to be 2.5 trillion 
US dollars, which is around 3-4% of the global 
GDP [3]. Some of the metal materials come in 
contact with acids used during _ prickling, 
descaling, and oil well acidizing [4, 5]. Corrosion 
is a natural process that cannot be eliminated but 
can be mitigated by adopting some approaches. 
These approaches include modification of 
metal/alloy, prevention of surface reactions 
(anodic and cathodic methods, protective coatings, 
modification of surface conditions, and the 
inhibition of surface corrosion (use of corrosion 
inhibitors) [6, 7]. The use of corrosion inhibitors is 
the most cost-effective and efficient way of 
mitigating corrosion. 

Corrosion inhibitor has been extensively used 
in mitigating corrosion in many industrial sectors. 
It is a liquid or solid substance which reduces the 
rate of corrosion when added to a suitable 
concentration without significantly changing the 
concentration of the corrosive agent [8, 9]. 
Corrosion inhibitors could be organic or inorganic 
inhibitors. The use of organic corrosion inhibitors 
is most preferred due to their high efficiency, less 
toxicity, and being eco-friendly and cheaply 
available [10, 11]. The efficacy of organic 
corrosion inhibitors usually depends on the 
presence of heteroatoms including oxygen, 
nitrogen, and phosphorous, and pi-electron which 
enables them to be adsorbed on the metal surface 
forming a barrier between the metal surface and 
environment and protecting the surface from an 
acidic environment [10, 12]. The adsorption of 
inhibitors on a metal surface can either be 
chemisorption (involving the interaction between 
the lone pairs of the molecule and the vacant d- 
orbital of the metal to form a coordinate bond) or 
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physisorption (involving the electrostatic 
interaction between the molecule and the 
electrically charged metal surface) [13]. 

Triazolopyrimidinone is a heterocyclic 
molecule that comprises triazole and fused 
pyrimidinone/4-hydoxylpyrimidine rings [14, 15]. 
Triazolopyrimidine compounds have been widely 
used in medicinal chemistry for their various 
therapeutic properties [16]. Triazolopyrimidine 
derivatives have been reported to possess 
antibacterial, antitumor, and anti-malaria 
activities. Due to the presence of heteroatoms (N, 
S that can donate an electron to the metal surface, 
fully conjugated pi-electron, its derivatives could 
act as potential anti-corrosive candidates. The 
corrosion inhibition efficiency of a compound can 
be effectively measured via experimental means. 
However, because of cost and time, theoretical 
methods are now employed to first investigate 
these molecules before laboratory experiments. 

Density functional theory (DFT) is a powerful 
tool employed to investigate the electronic 
properties of a molecule [17]. DFT provides 
insight into the transfer and distribution of 
electrons when molecules interact with a metal's 
surface. It has proven a very effective and accurate 
method of predicting the properties of molecules 
[18, 19], including corrosion inhibition potentials 
[20, 21]. Furthermore, the Monte Carlo (MC) 
simulation explains a molecule's orientation and 
mechanism of adsorption on the metal surface 
[22]. It predicts the adsorption energy when an 
inhibitor adsorbs on the surface of a metal. 
Therefore, this work investigated the corrosion 
inhibitive potentials of some triazolopyrimidinone 
derivatives (Fig. 1) using DFT and Monte Carlo 
simulation. These molecules were earlier reported 
by Istanbullu and coworkers as antimicrobial 
agents [16]. 


Methods 


Quantum chemical calculations 


Optimization of the structures of 
triazolopyrimidinone derivatives (A1-A6, B1-B6) 
was performed at the DFT/B3LYP/6-31G* level of 
theory on Spartan 14 software [23] having initially 
obtained the global minima (most stable 
conformers) [24]. All DFT calculations were done 
in the gas phase. The quantum parameters were 
calculated (equations 1-9) as described earlier 


[25]. 
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Fig. 1 Structures of the Triazolopyrimidinone derivatives under investigation. 


Al = 3H-5-(Chloromethyl)-2-(metylthio)-[1,2,4]triazolo[1,5-aļpyrimidine-7-one 

A2 = 3H-5-(Chloromethy])-2-(thiophene-2-yl)-[1,2,4]triazolo[1,5-a]pyrimidine-7-one 

A3 = 3H-5-(Chloromethy])-2-(4-chlorophenyl)-[1,2,4]triazolo[1,5-a]pyrimidine-7-one 

A4 = 3H-5-(Chloromethyl)-2-phenyl-[1,2,4]triazolo[ 1,5-a]pyrimidine-7-one 

A5 = 3H-5-(Chloromethy])-2-(4-nitrophenyl)-[1,2,4]triazolo[1,5-a]pyrimidine-7-one 

A6 = 3H-5-(Chloromethy])-2-(4-methylphenyl)-[1,2,4]triazolo[1,5-a]pyrimidine-7-one 

B1 = 3H-2-(Methylthio)-5-((4-(naphthalen-1-ylmethyl)piperazin-1-yl)methyl)-[1,2,4]triazolo[1,5 a]pyrimidin-7-one 
B2 = 3H-5-((4-(Naphthalen-1-ylmethyl)piperazin-1-yl)methyl)-2-(thiophen-2-yl)-[1,2,4]triazolo[1,5 a]pyrimidin-7-one 
B3 = 3H-2-(4-Chloropheny1)-5-((4-(naphthalen- 1-ylmethy])piperazin-1-yl)methyl)-[1,2,4]triazolo[1,5- a]pyrimidin-7-one 
B4 = 3A-5-((4-(Naphthalen-1-ylmethyl)piperazin-1-yl)methyl)-2-phenyl-[1,2,4]triazolo[1,5-a]pyrimidin-7-one 

B5 = 3H-5-((4-(Naphthalen-1-ylmethyl)piperazin-1-yl)methyl)-2-(4-nitrophenyl])-[1,2,4]triazolo[ 1,5-a]pyrimidin-7-one 
B6 = 3H-5-((4-(Naphthalen-1-ylmethyl)piperazin-1-yl)methyl)-2-(p-tolyl)-[1,2,4]triazolo[1,5-a]pyrimidin-7-one 
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molecules can be obtained using the Fukui 
function calculation. The Fukui function provides 
insight into the region where nucleophilic and 
electrophilic attacks occur in a molecule [26]. 
Charge of an atom in the neutral q(N), anionic 
q(N+1), and cationic g(N—1) states, that the Fukui 
functions can be uniquely calculated using 
equations 10-12. 


eq. 10 
fkt = [qk(N + 1)- qk(N)] (nucleophilic 
attack) eq. 11 
fi = [qk(N) — qk(N-1)] (electrophilic 
attack) eq. 12 


Afi(r) = fkt = f 


Monte Carlo simulation 


Compounds A1-6 and B1-6 were simulated using 
the Monte Carlo method to assess the adsorption 
and possible interactions between them and 
Fe(110) surface along the low-index miller facet 
[8, 22, 25]. Fe (110) plane was selected for its good 
stability and well-packed structure [27]. The 
molecules and Fe (110) were optimized using the 
Forcite module as earlier described [8, 25]. In a bid 
to mimic the experiment, the simulation was done 
in an aqueous environment by including water 
molecules. Adsorption of the molecules was done 
using the COMPASS force field in the adsorption 
locator module to obtain the adsorption energy 
(Eads) and the deformation energy (Eper) in 
kcal/mol [22]. 
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Results and Discussion 


Corrosion inhibition, frontier orbital energies, 
and reactivity descriptors 


The assessment of reactivity of the molecules is 
established on the frontier molecular orbitals 
energy (FMOs) theory of Fukui. The FMOs 
energies, that is, Emomo (energy of highest 
occupied molecular orbital) and ELumo (energy of 
unoccupied molecular orbital) are important 
quantum properties used to describe the reactivity 
of a molecular system. Eyomo and Epumo help in 
pointing out the orbitals that can donate and those 
that can accept electrons in a molecule. These 
energies help in predicting the corrosion inhibitory 
performance of molecules. Eyomo explains the 
ability of molecules to donate electrons. The 
HOMO energies ranged from (-6.42 to -5.63 eV) 
of the molecules follow the following trend: B6 > 
B4 = B2 > B1 > B5 > B3 > A6 > A2 > A4 > Al > 
A3 > A5 (Table 1). Molecule B6 has the highest 
HOMO energy (-5.63 eV) and the highest ability 
to donate electrons, that is the energy required by 
it to donate an electron is the lowest (ionization 
potential of 5.63 eV), while A5 has the lowest 
HOMO energy (-6.42 eV) and, hence, the highest 
number of ionization potential, making it the least 
ready to give up an electron [28]. The electron 
accepting ability of a molecule is described by the 
EruĮmo value, low Erumo value favors a high 
potential of electron acceptance [7]. The ELumo of 
the molecules follow the following trend: B1 > A1 
> B2 > B6 > B4 > A2 > A6 > A4 > B3 > A3 > B5 
> A5. Molecule B1 has the highest LUMO (-0.94 
eV), hence the lowest ability to accept an electron 
that is its electron affinity is the lowest (0.94 eV), 
while A5 has the lowest LUMO (-3.04 eV), hence 
the highest electron affinity (3.04 eV) [28]. The 
Energy band gap (Ev) explains the reactivity and 
stability of a molecule during adsorption on a 
metal surface. A lower value in Eg, indicates an 
increase in the reactivity of a molecule, enabling 
improved inhibitory effects while the higher value 
of E indicates high chemical stability, leading to 
a decrease in inhibitory effect [29]. The Eg values 
for the studied compound are of the order: Al > 
B1 > A2 > A6 > A4 > A3 > B2 > B6 > B4 > B3 > 
AS > B5 (Table 1). Generally, the molecular 
orbitals of molecules B1-B6 are more stabilized 
over their A series due to the higher molecular size 
of each counterpart. This also led to a decrease in 
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Table 1 Electronic parameters and reactivity descriptors of the compounds: frontier molecular orbital energies (E nomo and 
ELumo), ionization potential (I), electron affinity (A), electronegativity (y), chemical hardness (nV), chemical softness (ò), 
electrophilicity (œ), fraction of electron transferred (AN) and electron back-donation (AEback-donation). All values were 


calculated at the B3LYP/6-31G(d) level of theory. 


Exomo ELumo Eg I A ya n a ji 
Compounds (eV) (eV) (eV) (eV) (eV) (eV) (eV) (eV"!) (eV) AN AEback-donation 
Al -6.20 -1.30 4.90 6.20 1.30 3.75 245 0.41 2.87 0.66 -0.61 
A2 -6.12 -1.70 4.42 6.12 1.70 3.91 2.21 0.45 3.45 0.69 -0.55 
A3 -6.25 -2.04 4.21 625 204 4.15 2.11 0.48 4.08 0.68 -0.53 
A4 -6.13 -1.82 4.31 6.13 1.82 3.97 2.15 0.46 3.66 0.70 -0.53 
AS -6.42 -3.04 3.38 642 3.04 4.73 1.69 0.59 6.61 0.67 -0.42 
A6 -6.08 -1.72 4.36 6.08 1.72 3.90 2.18 0.45 3.48 0.71 -0.54 
B1 -5.66 -0.94 4.72 5.66 0.94 3.30 2.36 0.42 2.30 0.78 -0.59 
B2 -5.64 -1.49 4.15 5.64 149 3.56 2.07 0.48 3.06 0.82 -0.51 
B3 -5.68 -1.86 3.82 5.68 1.86 3.77 1.91 0.52 3.70 0.84 -0.47 
B4 -5.64 -1.60 4.04 5.64 1.60 3.62 2.02 0.49 3.24 0.83 -0.50 
B5 -5.71 -2.94 2.77 5.71 2.94 4.32 1.38 0.72 6.74 0.96 -0.34 
B6 -5.63 -1.51 4.12 5.63 1.51 3.57 2.06 0.49 3.09 0.83 -0.52 


the energy band gaps of molecules B1-B6 over 
their corresponding A analogs. For example, A5 
has 3.38 eV while its analog B5 has 2.77 eV. From 
both series, the nitro group resulted in the lowering 
of the energy band gap, as was observed in the 
literature [30], but was lower in A5 as a result of 
higher molecular size. Compound B5 has the 
lowest band gap energy indicating that it is the 
most reactive, therefore could possess better 
corrosion inhibition potentials than other 
compounds under investigation. The low energy 
band gap values of A5 and B5 are due to the 
attachment of electron-withdrawing nitro groups 
in the structures; nitro-substituted molecules are 
usually characterized by low energy band gaps 
[30]. 

Global reactivity descriptors are important 
parameters used in describing and predicting the 
inhibitory effect of potential inhibitors. The 
electronic properties and global reactivity 
descriptors of the molecules are presented in Table 
1. A maximum ionization potential signifies that 
the compound has a low electron donating ability, 
while a low ionization potential signifies a high 
electron acceptance ability [12, 31]. Compounds 
B6, B4, and B2 displayed the lowest ionization 
energy (5.63, 5.64 and 5.64 eV, respectively) 
indicating the high electron-releasing potential of 
the compounds. Electronegativity (%) is the ability 
of an atom in a molecule to attract electrons toward 
itself [32]. Lower electronegativity values enable 
molecules to easily accomplish electron 
equilibrium, increasing their reactivity. The 


electronegativity values of the studied compounds, 
as shown in Table 1, are in the range of 3.30 eV — 
4.32 eV with compound G showing the lowest 
electronegativity value. The global hardness (n) 
and softness (ò) are determined in compliance with 
the HSAB theory [33]. The theory says that the 
compounds prefer to bind with systems where they 
pose the same or similar hardness (softness) and 
the rise in interaction energy occurs with a 
decrease in hardness and increase in softness [34]. 
Since hardness is directly from the energy band 
gap, B3, A5, and B5 have the lowest chemical 
hardness and chemical softness. (1.91 eV, 1.69 eV 
and 1.38 eV) and softness value (0.52 eV, 0.59 
eV! and 0.72 eV). 

Energy stabilization via the flow of electrons 
is measured by the value of global electrophilicity. 
Nucleophiles possess low _ electrophilicity. 
Molecule B1 is the best nucleophile while B5 is 
the best electrophile. Furthermore, molecules are 
absorbed into the metal surface by donating 
/transferring electrons. Transferring more 
electrons enhances inhibition efficiency. The 
electron transfer (AN) number indicates better 
interaction between the molecule and metal 
surface since they are all less than 3.6 [25]. The 
AN values shown in Table 1 follow the order: B5 
> B3 > B4 = B6 > B2 > B1 > A6 > A4 > A2 > A3 
> A5 > Al. Negative values of AEback-donation 
displayed by all molecules indicate that charge 
transfer between the surface of the metals and the 
molecules is energetically favored, implying that 
all molecules can make electrons available via 
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Fig. 2 HOMO map (A), LUMO map (B) and electrostatic potential map (C) of triazolopyrimidinone derivative A1. 


back-donation [35]. 


Surfaces (HOMO, LUMO and electrostatic 
potential maps) 


The HOMO, LUMO and electrostatic potential 
(ESP) maps of compound A1 are displayed in Fig. 
2. Figs. S1-S11 show the HOMO, LUMO and 
electrostatic potential (ESP) maps of compounds 
A2, A3, A4, A5, A6, B1, B2, B3, B4, B5 and B6. 
The HOMO of Al spreads across the entire 
molecule, leaving out the chloromethyl and the 
methyl on the thiomethyl groups (Fig. 2a), while 
the LUMO covers the same regions, adding the 
chloromethyl group and leaving out the entire 
thiomethyl group (Fig. 2b). The ESP map of A1 
indicates that the carbonyl oxygen (red) on the 
pyrimide-7-one and chloromethyl (yellow and 
orange) regions are the sites that will mostly favor 
electrophilic attack while the triazole hydrogen, 
sulfur on methylthiol group, triazole and 
pyrimidine ring favor nucleophilic attack (Fig. 2c). 
The carbonyl oxygen and the triazole hydrogen are 
the best sites for electrophilic and nucleophilic 
attacks, respectively. The HOMO map of 
compounds A2 covers the triazole and pyrimidine 
rings (Fig. Sla), while the LUMO covers the 
thiophene-2-yl and parts of the triazole and 
pyrimidine rings, leaving out the chloromethyl 
group (Fig. S1b). The ESP map of A2 indicates 
that the carbonyl oxygen (red) on the pyrimide-7- 
one and chloromethyl (yellow and orange) regions 
are the sites that will mostly favor electrophilic 
attack while the triazole hydrogen, entire 
thiophene group, triazole and pyrimidine ring 
favor nucleophilic attack (Fig. S1c). The HOMO 
map of compound A3 covers the entire molecule 
except for the chlorophenyl and chloromethyl 
groups (Fig. S2a) while the LUMO covers majorly 
the chlorophenyl and triazole rings (Fig. S2b). The 
ESP map of A3 indicates that the carbonyl oxygen 


(red) on the pyrimide-7-one is the site that will 
mostly favor electrophilic attack while other 
regions will favor nucleophilic attack (Fig. S2c). 
The HOMO map of compounds A4 covers the 
triazole and pyrimidine rings (Fig. S3a), while the 
LUMO covers the phenyl, triazole and part of the 
pyrimidine ring (Fig. S3b). The ESP map of A4 
indicates that the carbonyl oxygen (red) on the 
pyrimide-7-one and the chloromethyl (orange and 
yellow) are the sites that will mostly favor 
electrophilic attack while other regions will favor 
nucleophilic attack, with the most preferred being 
the triazole hydrogen (Fig. S3c). The HOMO map 
of compound A5 covers the triazole and 
pyrimidine rings (Fig. S4a), while the LUMO 
covers the nitrophenyl ring (Fig. S4b). The ESP 
map of A5 indicates that the carbonyl oxygen (red) 
on the pyrimide-7-one, nitro group on the 
nitrophenyl ring (yellow) and the chloromethyl 
(orange and yellow) are the sites that will mostly 
favor electrophilic attack while other regions will 
favor nucleophilic attack, with the most preferred 
being the triazole hydrogen (Fig. S4c). The HOMO 
map of compound A6 covers the triazole and 
pyrimidine rings (Fig. S5a), while the LUMO 
covers the entire molecule (Fig. S5b). The ESP 
map of A6 indicates that the carbonyl oxygen (red) 
on the pyrimide-7-one and the chloromethyl 
(orange and yellow) are the sites that will mostly 
favor electrophilic attack while other regions will 
favor nucleophilic attack, with the most preferred 
being the triazole hydrogen (Fig. S5c). 

For B1-B6, the HOMO and LUMO of BI 
spread across the naphthalene-1l-ylmethyl group 
(Figs. S6a and b). The ESP map of B1 indicates 
that the triazole and pyrimidine rings (red) are the 
sites mostly preferred for the electrophilic attack 
while the triazole hydrogen and other regions favor 
nucleophilic attack (Fig. S6c). The HOMO of B2 
covers the naphthalene-1-ylmethyl group (Fig. 
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Table 2 Charges on the anionic, neutral and cationic species qk (N+1), qk (N) and qk (N-1), respectively. The Fukui 
charges for nucleophilic (fk*) and electrophilic attack (fk) and the Fukui function (Afx) on triazolopyrimidinone derivative 


Al. 

Atom qx (N+1) qk qx (N-1) fkt fk Afk 
Cl 0.366 0.326 0.315 0.04 0.011 0.029 
C2 -0.305 -0.276 -0.193 -0.029 -0.083 0.054 
C3 0.508 0.622 0.653 -0.114 -0.031 -0.083 
Ol -0.562 -0.494 -0.345 -0.068 -0.149 0.081 
NI -0.451 -0.442 -0.405 -0.009 -0.037 0.028 
C4 0.787 0.818 0.845 -0.031 -0.027 -0.004 
N2 -0.598 -0.56 -0.469 -0.038 -0.091 0.053 
C5 -0.332 -0.436 -0.446 0.104 0.01 0.094 
Cll -0.729 -0.075 0.01 -0.654 -0.085 -0.569 
N3 -0.289 -0.275 -0.202 -0.014 -0.073 0.059 
C6 0.32 0.329 0.332 -0.009 -0.003 -0.006 
N4 -0.714 -0.708 -0.672 -0.006 -0.036 0.03 
Sl 0.126 0.181 0.393 -0.055 -0.212 0.157 
C7 -0.577 -0.588 -0.66 0.011 0.072 -0.061 


S7a) while the LUMO covers the thiophene, 
triazole and part of the pyrimidine rings (Fig. S7b). 
The ESP map of B2 indicates that the major 
electrophilic sites are the carbonyl oxygen (red) on 
the pyrimide-7-one with other parts in yellow and 
orange also favors electrophilic attack while other 
parts favor nucleophilic attack (Fig. S7c). The 
HOMO of B3 covers the same region as in B2 (Fig. 
S8a), while the LUMO covers the chlorophenyl, 
triazole and some parts of the pyrimidine rings 
(Fig. S8b). The ESP map of B3 indicates that the 
carbonyl oxygen (red) on the pyrimide-7-one is the 
site that will mostly favor electrophilic attack 
while other regions will favor nucleophilic attack 
(Fig. S8c). The HOMO of B4 covers the same 
region as in B2 and B3 (Fig. S9a), while the 
LUMO covers the phenyl, triazole and part of the 
pyrimidine ring (Fig. S9b). The ESP map of B4 
indicates that the carbonyl oxygen (red) on the 
pyrimide-7-one and part of the pyrimidine ring 
(yellow) are the sites that will mostly favor 
electrophilic attack while other regions (the 
phenyl, triazole and naphthalene-1-ylmethyl 
groups) favor nucleophilic attack, with the most 
preferred being the triazole hydrogen (Fig. S9c). 
The HOMO of B5 covers some part of the 
naphthalene-1-ylmethyl group (Fig. $10a), while 
the LUMO covers the nitropheny] ring (Fig. S10b). 
The ESP map of B5 indicates that the carbonyl 
oxygen (red) on the pyrimide-7-one, some parts of 
the pyrimidine ring and the naphthalene-1- 
ylmethyl group (yellow, red and orange) favor 
electrophilic attack while other regions will favor 


nucleophilic attack (Fig. S10c). The HOMO of B6 
covers the naphthalene-l-ylmethyl group (Fig. 
Sila), while the LUMO covers the 
methylbenzene, triazole and some parts of the 
pyrimidine rings (Fig. S11b). The ESP map of B6 
indicates that the carbonyl oxygen (red) on the 
pyrimide-7-one and some parts of the naphthalene 
ring (yellow) are the sites that will mostly favor 
electrophilic attack while other regions will favor 
nucleophilic attack, with the most preferred being 
the triazole hydrogen (Fig. 11S5c). The molecules 
all have sites that can accept and donate electrons, 
implying that they can bind to the surface of the 
metal via donation and back-donation of electrons, 
a rationalization of their probable use as corrosion 
inhibitors. 


Local reactivity descriptors and 


distribution 


charge 


The reactive sites on the compounds could be 
determined by the Mulliken charge distribution. 
Equations 10-12 show the calculated charges at the 
cationic, neutral and anionic levels (Tables 2 and 
S1—S11). Electrophilic atoms are characterized by 
positive Af, values, while nucleophilic atoms are 
characterized by negative Af, values. The favorite 
nucleophilic site for molecule Al is atom C5 
(Table 2) with the highest value of fr” (0.104) and 
a positive Af, value (0.094), while the favorite 
electrophilic site is atom C7 with fw (0.072) and a 
negative Af, value (-0.061). The favorite 
nucleophilic site for molecule A2 is atom C5 
(Table S1) with the highest value of fx* (0.104) and 
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Fig. 3 Side view (A) and top view (B) of adsorption of triazolopyrimidinone derivative Alon Fe(110). 


a positive Af, value (0.096) while the favorite 
electrophilic site is atom C7 with fw (0.023) and a 
negative Af, value (-0.013). The favorite 
nucleophilic site for molecule A3 is atom C5 
(Table S2) with the highest value of fx* (0.103) and 
a positive Af, value (0.096), while the favorite 
electrophilic site is atom C7 with fw (0.01) and a 
negative Af, value (-0.008). The favorite 
nucleophilic site for molecule A4 is atom C5 
(Table S3) with the highest value of fx* (0.104) and 
a positive Af, value (0.095), while the favorite 
electrophilic site is atom C7 with fx (0.07) and a 
negative Af, value (-0.006). The favorite 
nucleophilic site for molecule A5 is atom C8 
(Table S4) with the highest value of fr”? (0.447) and 
a positive Af, value (0.46) while the favorite 
electrophilic site is atom C7 with fx (0.443) and a 
negative Afk value (-1). The favorite nucleophilic 
site for molecule A6 is atom C5 (Table S5) with 
the highest value of fk* (0.105) and a positive Afk 
value (0.99) while the favorite electrophilic site is 
atom C13 with fx (0.009) and a negative Afk value 
(-0.006). 

For B1-B6, the favorite nucleophilic sites for 
molecule B1 are atoms N1 and N4 (Table S6) with 
the highest values of fr” (0.022 each) and positive 
Afk values (0.024 and 0.025, respectively), while 
the favorite electrophilic site is atom C5 with fx 
(0.074) and a negative Af, value (-0.069). The 
favorite nucleophilic site for molecule B2 is atom 
NI (Table S7) with the highest values of fkt 
(0.015) and a positive Af, value (0.026), while the 
favorite electrophilic sites are atoms C5 and C21 
with fx value of (0.025 each) and negative Afk 


values (-0.022 and -0.025, respectively). The 
favorite nucleophilic site for molecule B3 is atom 
N1 (Table S8) with the highest value of fx* (0.02) 
and a positive Af, value (0.02) while the favorite 
electrophilic site is atom C10 with f, value of 
(0.539) and a negative Af, value (-0.575). The 
favorite nucleophilic site for molecule B4 is atom 
NI (Table S9) with the highest value of fx* (0.018) 
and a positive Af, value (0.02), while the favorite 
electrophilic site is atom C15 with fw value of 
(0.369) and a negative Af, value (-0.771). The 
favorite nucleophilic site for molecule B5 is atom 
C9 (Table S10) with the highest value of fi? 
(0.067) and a positive Afk value (0.111), while the 
favorite electrophilic site is atom C5 with fx value 
of (0.067) and a negative Af, value (-0.062). The 
favorite nucleophilic site for molecule B6 is atom 
NI (Table S11) with the highest value of fi? 
(0.017) and a positive Afk value (0.02), while the 
favorite electrophilic site is atom C5 with fx value 
of (0.483) and a negative Af, value (-1.036). In 
essence, the nucleophilic and electrophilic sites on 
the corrosion inhibitors are located on different 
carbon atoms for effective metal-inhibitor 
interactions. 


Monte Carlo simulation 


The total energy of iron-inhibitor (Erot) and 
adsorption energy (Eads), which is a sum of the 
rigid adsorption energy (Erigia) and the 
deformation energy (Eper) were calculated (Table 
3). Erot, a measure of stability of the substrate- 
adsorbate systems is in the order: Fe(110)-A6 < 
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Table 3 Monte Carlo simulation results for Fe(110)-inhibitor systems for all compounds (All energies are in kcal/mol). 
Eror = Total energy, Eaas = Adsorption energy, Erig = Rigid adsorption energy, Epef = Deformation energy. 


System Ertot Erig Eper 
Fe(110)-A1 -11266.02 -13257.7 -11553.4 -1704.3 
Fe(110)-A2 -11276.45 -13281.3 -11584.8 -1696.5 
Fe(110)-A3 -11290.40 -13290.1 -11604.3 -1685.7 
Fe(110)-A4 -11188.26 -13174.6 -11460.3 -1714.2 
Fe(110)-A5 -11106.14 -13095.0 -11393.7 -1701.2 
Fe(110)-A6 -11337.98 -13313.2 -11600.0 -1713.2 
Fe(110)-B1 -11323.03 -13231.6 -11519.3 -1712.2 
Fe(110)-B2 -11101.56 -13026.4 -11299.1 -1727.3 
Fe(110)-B3 -11143.26 -13045.3 -11332.9 -1712.3 
Fe(110)-B4 -11247.01 -13166.5 -11451.4 -1712.2 
Fe(110)-B5 -11331.17 -13235.1 -11527.6 -1707.4 
Fe(110)-B6 -11243.56 -13142.7 -11422.8 -1719.8 

Fe(110)-B5 < Fe(110)-B! < _ Fe(110)-A3 < and the metal surface via back-donation. The Etot 


Fe(110)-A2 < Fe(110)-A1 < Fe(110)-B4 < 
Fe(110)-B6 < Fe(110)-A4 < Fe(110)-B3 < 
Fe(110)-A5 < Fe(110)-B2. The values displayed 
for Ett show that the adsorption of molecules to 
the metal surface is thermodynamically favored. 
Low adsorption Eaas, as observed in all 
compounds is associated with a strong adsorption 
of the inhibitor molecules on Fe(110) surface [36, 
37]. Fig. 3 and S12-22 show the top and side views 
of the adsorbed molecules on the surface of 
Fe(110) in an aqueous solution, with all molecules 
arranged parallel on the plane of Fe(110) and 
providing good coverage on the surface. This 
shows that the molecules can bind to the metal 
surface via covalent and non-bonded interactions 
[38], buttressed by the negative values of AEback- 
donation [28]. 


Conclusions 


The corrosion inhibition potentials of some 
derivatives of triazolopyrimidinone were 
investigated via the density functional theory 
approach and Monte Carlo simulations. The 
electronic properties and global and local 
reactivity descriptors were calculated to ascertain 
the reactivity, stability and reactive sites of the 
molecules. Adsorption locator calculation was 
done to obtain the total energy and adsorption 
energy of the molecules on Fe(110) surface. The 
results revealed the existence of donor and 
acceptor parts in the molecules. Low energy band 
gap values were observed for the nitro substituted 
derivatives, as observed in the results of previous 
researchers. Negative AEpack-donation Values in the 
molecules suggest charge transfer between them 


values suggest thermodynamically favored 
adsorption between the molecules and the metal 
surface. These molecules could be used as anti- 
corrosive agents. 
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